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Recently Flach et.al.0| used a symmetry analysis to predict the appearance of directed energy 
current in homogeneously spatially extended systems coupled to a heat bath in the presence of an 
external ac field E{t). Their symmetry analysis allowed them to make the right choice of E{t) so 
as to obtain symmetry breaking which causes directed energy transport for systems with a nonzero 
topological charge. Their numerical simulations verified the existence of the directed energy current. 
They argued that the origin of their strong rectification in the underdamped limit is due to the 
excitation of internal modes and their interaction with the translational kink motion. The internal 
mode mechanism as a cause of current rectification was also proposed by Salerno et.al.Q. We use 
a rigorous collective variable for nonlinear Klein-Gordon equations to prove the rectification of the 
current is due to the excitation of an internal mode, V{t), which describes the oscillation of the slope 
of the kink, and to a dressing of the bare kink by the ac driver. The internal mode T{t) is excited 
by its interaction with the center of mass of the kink, X{t), which is accelerated by E{t). The 
external field E{t) also causes the kink to be dressed. We derive the expressions for the dressing and 
numerically solve the equations of motion for T{t), X(t), and the momentum, P{t), which enable 
us to obtain the explicit expressions for the directed energy current and the ac driven kink profile. 
We then show the directed energy current vanishes unless the slope T{t) is a dynamical variable and 
the kink is dressed by the ac driver. 

PACS numbers: 05.45.-a, 05.45.Yv 



I. INTRODUCTION 

In a recent paper S. Flach et.al.l'l studied the appear- 
ance of directed energy currents in homogeneous spatially 
extended systems described by nonlinear field equations 
coupled to a heat bath in the presence of an external 
ac field E{t). As pointed out in rectifying energy 
transform using fluctuations has been studied in connec- 
tion with such problems as molecular motors in biological 
systems |3j] , electrical currents in superlattices A, A A Zl , 
voltages in Josephson junction coupled systems |a.l9l[l0l| 
and other problems. 

The authors of Ref.[l] showed by a symmetry analy- 
sis that the correct choices of E{t) lead to directed en- 
ergy transport for nonlinear Klein-Gordon systems with 
a nonzero topological charge. They used numerical sim- 
ulations of the ac driven Klein-Gordon equation which 
confirmed their predictions which generalized recent rig- 
orous theories of currents generated by broken time-space 
symmetries to the case of interacting many particle sys- 
tems ^3 . They did this by replacing the fluctuations 
as a superposition of ac driving fields and uncorrelated 
white noise. They also showed the persistence of directed 
currents in the Hamiltonian limit of systems exposed to 
ac fields but decoupled from the heat bath. The authors 
of Ref.[l] then argued that the origin of the observed 
strong rectification in the underdamped limit is due to 
the nonadiabatic excitation of internal kink modes and 
their interaction with the translational kink motion. 

In this paper we use a rigorous collective variable (CV) 
theory for nonlinear Klein Gordon equations derived in 
Refs.[13,14] to prove that an external ac field causes the 
CV's for the center of mass X(t) and the slope T{t) to 
become time dependent and to interact with each other. 



The ac driver in addition to inducing time dependence 
in X{t) the center of mass of the kink and the slope V{t) 
causes the kink to be dressed by phonons. The dressing 
changes the shape of the kink and increases the coupling 
of r(t) to X{t). We show the non- vanishing of the dress- 
ing is a necessary condition for breaking time inversion 
symmetry of the energy current. However, the dressing 
of the kink alone in the absence of a time dependent col- 
lective variable, F(t), cannot cause current rectification. 
Consequently, we prove that the existence of the time de- 
pendence of r(i) and the dressing are necessary for 
time inversion symmetry breaking. 

In Sec. II we derive the equations of motion for X{t) 
and F(i) including the terms due to the dressing of the 
kink by phonons. We present our results for the solutions 
X{t) and T{t) and for the generation of directed energy 
currents in Sec. Ill and in Sec. IV and we discuss our re- 
sults. The derivation of the dressing x is given in the 
appendix. 



II. DERIVATION OF CV EQUATIONS OF 
MOTION 

Before deriving the CV equations of motion used in this 
paper, we will make a few remarks about CV treatments 
of the Klein-Gordon equations. The first approach which 
is derived in Refs.[13,14] and used in this paper, is to treat 
the center of mass X{t) and the slope T{t) as collective 
variables which satisfy coupled second order differential 
equations, which also depend on the dressing of the kink. 
In this approach the equations of motion for X{t) and 
T{t) are not manifestly relativistic invariant. What has 
been proven is that when the solutions X(t) and T{t) are 
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inserted in the kink (j)[X{t),T{t)] that (f> satisfies the rel- 
ativistic invariant nonUnear Klein-Gordon equation. An 
analogous well-known example of a non-manifestly rel- 
ativistic case is the use of the non-relativistic Coulomb 
gauge which leads to the relativistic solution E{t) and 
B{t) of Maxwell's equations. The second approach would 
be to consider a single CV, X{t). The equation of motion 
for X{t) is fourth order in time because the Lagrangian 
contains the second derivative X{t). The Lagrangian of 
a theory that contains a second derivative, X{t), leads 
to equations of motion that contain the fourth deriva- 
tive 'X{t). The fourth order equation for X{t) has the 
same number of degrees of freedom as the equivalent CV 
theory for two variables which consists of two coupled 
second order equations for X{t) and r{t). 

We outline the derivation of the equations of motion 
for the collective variables X{t) and r{t) which are de- 
rived in detail in Refs.[13,14]. The damped nonlinear SG 
equation for the field in the presence of an external 
potential, V{(j)), is 



- sin (j) + /3(f>j, 



dV 



(1) 



where f3(f>^t is the damping due to the heat bath, and 
where we are using dimensionless variables where the ve- 
locity of the phonons is c = 1 . We introduce the collective 
variables by writing the solution cf) of Eq.(l) in the form 



(2) 



where ^ = r{t)[x — X{t)] and the single kink solution 
am] is 



a[^{t)] ^ 4:tan-^exp[r{t){x - X{t))] 



(3) 



and x[^(i)] is the dressing of the kink by phonons due 
to the external potential ¥{(/)) which for the applied ac 
field of this paper is given by 



V{(j)) — ( eiCOSLut+e2COs[2Lut+9 



(4) 



9 is an arbitrary phase and ei, and £2 are perturbation 
parameters i.e., we solve for to first order in ei and 
£2- The CV's are the center of mass X{t) and the slope 
of the kink evaluated at its center is 2r(t). In Ref. [2] a 
directed kink motion for the SG was obtained numerically 
for the first time using V{4i) in Eq.(4) for a wide range of 
momenta with an analytic approach for small momenta. 

The equations of motion for X and T each contain 
many terms proportional to integrals of x, its time deriva- 
tives and spatial derivatives, x is a solution of the lin- 
earized ac driven SG equation and is proportional to £1 
and £2. In the appendix we solve for x- The solution for 
X is 



where 



sech^^ 



(5) 



/(<) = (£i/2)costjt 



1-LU 



l+UJ 



+ (£2/2) cos(2Ljt -t- 6*) 



/32 + (l-Cj)2 /32 + (l+t^)2 

1 - 2w 



/32 + (1 _ 2cj)2 
1 + 2lu 

+ {1 + 2ujy 



mi 



Since x is an even function of ^ many of the terms in 
Eqs.(2.4a) and (2.5a) of Ref.[13,14] for X and f that 
depend on integrals of x vanish. The only terms which 
survive are 



(1 - bx)M 



X 



X + X{t/T) + l3X 



2^/1 

+ T^<j'\x"){l-X') 

- it/T)^a'\ex") 

2(f/r)(///)(a'|^x) 

- (f/r)(a'|ex) (6) 

where Mx = T{a'\a') = 8T , bx = {r/Mx){a"\x) = 
and where {f\g) = J The corresponding 

equation for F is 



J 



{l-br)Mr 



f - 3rV2r -f (Mx/2r) (1 ~x^)+(3t 



= i2Xr/r'){(a'\x''0 + 2iX/m/ma'\x') + {X/mcT'\x') (7) 

I 



where Mr = T-^{^<j'\S,a') = (27r2/3r3) and br = 
(r3/Mr)-i(eV'|x) = 0. 

We next eliminate the F term in Eq. (6) and the X term 
in Eq.(7) by using the zeroth order in £1, and €2 equations 
for X and F. The elimination is justified because the 
corresponding terms are both multiplied by x which is 
already first order in €1, and £2. 



The zeroth order expression for X is X ^ — X(f/F) 
and for f is f = 3f 2/2F - Afx(2rMr)"i(l - X^). When 
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we substitute for X in Eq.(7), we obtain 
f + (3t = (3f V2r) + (6/7r2)r3(l - 

+ (3/27r2) 8rx{2{^a'\x") 

and when we substitute for F in Eq.(6) we obtain 



(8) 



X + X{t/T) +I3X = 27r/i - 2/(f /r)(a'^|x') 
+T\l-X')(^{a'\x")-i6/7r'){a'\^x') 

-(f/r)2f(a'|^\')-3/2(<7'C|x')) (9) 



The momentum P conjugate to X is P = MxX = STX. 
Consequently we can write Eq.(8) for X in terms of P 
i.e, 



dP 
dt 



+ /3P = 27r/i-2/(r/r)(a'^|x'> 



+(r^-^)((-'lx")-;^(-'l^x')) 



-C-f[{a'\e^')-m^'iW) 



and Eq.(8) for f in terms of P becomes 



(10) 



f - (3f V2r) - (6/7r2)r 



1 - + (p/8)^ 



+ /3r 



(3/27r2) 



Pt{2{ia'\x")-{ia'\x!) 



+TP{f/ma'\x!\ 



(11) 



Finally after evaluating the integrals and replacing X"^ 
by we obtain the final form of our equations of 

motion for P{t) and T{t) 



P + l3P-2nh = fit) 



0.47(8/7r)(f/r)2 

(8/3) (f/r) 

(8/7r2)fr2-(P/8)2 



(12) 



and 



r + /?r - 3r72r - (6/77^ r 1 - + (p/s) 



,2\-l 



5/2-7rVl6 Pf- /(t)//(t) Pr 



(13) 



Before solving the equations of motion for P and T it 
is worth making a few remarks about the properties of 
the coupled equations. In Eq.(12) the ac driver, /i(i), 
directly drives P while in Eq.(13) for T the ac driver, 
fi{t), does not directly drive T because the F mode,|f , 
is orthogonal to the ac driver. However, both P and F 
see the ac driver indirectly through the dressing x which 
is proportional to f{t), Eq.(A.7) which also depends on 
the two frequencies cj and 2a;. As long as w < 0.5 the 
phonon radiation is small. In this paper we consider only 
frequencies which are much less than one, which is the 
beginning of the lower band edge in the miits of this pa- 
per. Consequently in this paper the emission of phonons 
is negligible. Since there are no modes of X{t) and T{t) 
in the band gap of the SG, there is no excitation of inter- 
nal gap modes by the ac driver as there is e.g., in (j)^ and 
the double sine-Gordon. The dressing x which is propor- 
tional to e changes the shape modes to (t.x + x.x and to 
c,r + X,r in addition to changing the frequency of P and 
F directly. 

Before discussing the results we give a discussion of 
the symmetries of the coupled Eqs.(12) and (13). The 
first symmetry is rcfcred to as the shift symmetry of 
the driver which is P ^ — P and t ^ t + ^ provided 
hit) = -hit + §) and fit) = -fit + §) are always 
shift symmetric if and only if a Fourier expansion con- 
tains only odd terms. Thus /i(t) and fit) in Eqs.(12) 
and (13) always violate shift symmetry. A second sym- 
metry is time inversion symmetry i.e., P — P when 
t ^ —t and /? = 0. Eqs.(12) and (13) satisfy time inver- 
sion symmetry when /3 = and ^ = 0, ±mr. When /3 is 
small, time inversion is approximately satisfied. 



III. RESULTS OF SIMULATIONS 

In this section we present the computer solutions of 
Eqs.(12) and (13) for F(f), Pit), and for the time average 
of the energy current. In our units the energy current 
Jit) is equal to P(f) because 



At) 



'^AO 'jAC)dx = srx = Pit) 



where a[^] =4 tan ^ exp[ T{t) (a;— X(t)) ] . Consequently 
the time average of J(t), (Jit)), is equivalent to the time 
average of P(t), (P(t)). 

In Fig.(l) we show the results for (Pit)) for a range of 
values LU, /3, ei, and e2 which show clearly the directed 
energy current as a function of 9. We see that there is 
symmetry breaking for all the sets of parameter values. 

In Fig. (2) we show (P(t)) as a function of 9 for fixed 
values of u> and ei = e2 and various values of /3. If /3 ^ 
then time inversion symmetry is not valid. However as /3 
goes to zero, time inversion symmetry is approximately 
restored at ^ = 0, ±mr. Thus we observe for small (3 as 
(3 decreases, exactly the same behaviour as in Fig.(l) of 
Ref. [1] that is the smaller /3 the larger the value of (P(0 = 
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FIG. 1: This figure demonstrates symmetry breaking,i.e., the 
nonvanishing of (P{t)} as a function of 9 for various values of 
the parameters a;, e and (3. Solid curve: w = 0.1, ei = €2 = 
0.03 and /3 = { The curve is multiplied by 0.15), dashed 
curve: to = 0.3, ei = 0.3, €2 ~ ei/\/3 and /3 = 0.2, dotted 
curve: cu = 0.1, ei = £2 = 0.05, and = 0.12, dash-dotted 
curve: w = 0.25, ei = 0.16, £2 = ei/V2 and /J = 0.15. {P{t)) 
has the units of momentum and 6 is in radians. 
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FIG. 2: {P{t)) as a function of 6* for a; = 0.1, ei = €2 = 
0.03 for various values of /3 show a monotonic decrease of the 
amplitude of {P{t)) as the damping /3 increases. Solid curve: 
P = 0.02, dashed curve: [3 = 0.05, dotted curve: (3 = 0.12. In 
this simulation {P{t)} decreases as /3 increases but the values 
of {P{t)) are so small that they can not be distinguished on 
the figure. {P{t)} has the units of momentum and 6 is in 
radians. 



0)), but the smaller the value of at which {P{0)) = 0. 
The values in Fig.(2) at 6I = are (P(i)) = 0.0008 when 
(3 = 0.12, {P{t)) = 0.002 when /? = 0.05 and (P(i)) = 
0.005 when /3 = 0.02. 

In the three typical examples of P{t) and T{t), Figs. (3), 
(4) and (5) we have selected the shapes and magnitudes 
vary considerably. The shapes of the P{t) and T{t) curves 
show the effect of being driven by an ac driver with two 
frequencies which causes the curves to vary in amplitude 
and shape when we vary the parameters to, ei, 62, and /3. 
The changes in shape of T(t) arc often striking because 
r(i) is a very nonlinear oscillator which has a complicated 
response to the ac driver and the dressing %. Whereas 
in lowest order, the equation for P(t) is linear. Gener- 
ally the magnitudes of both T{t) and P{t) increase with 



increases in the strength of ei, and 62. The variable T{t) 
oscillates about an average value of (r(t)) which is greater 
than Fq = 1, the unperturbed kink value of F. It usually 
also takes instantaneous values less than one. The rel- 
ative change in slope, AF/F, varies from a few percent 
to as much as one hundred percent and is strongly de- 
pendent on the magnitude of ei and £2- Large values of 
AF /F represent large distortions of the shape of the kink. 
When we compare P{t) and F(i) for the kink dressed by 
X with the bare kink we find a strong dependence on the 
phase which leads to different shapes and amplitudes of 
P{t) and r{t) for different 9. For example in Fig. (5) for 
u, = 0.1, (3 = 0.02, ei = £2 = 0.03 and 9 near 9 = 0, ±mr 
the slope F(t) of the dressed kink is a pattern of single 
peaks while for the bare kink with the same parameters 
T{t) is a pattern of double kinks. On the other hand for 
9 appreciably different from 6* = 0, ±n7r e.g., 9 = 1.61 — tt 
the differences in shape of P{t) between the bare and 
dressed kinks are relatively minor. At the same time al- 
though the shapes of F(i) arc qualitatively similar the 
bare kinks have appreciably reduced amplitudes of T{t). 

IV. DISCUSSION 

The ac driver causes the center of mass X and the 
slope F to become time dependent and the kink to 
be dressed by phonons given by the expression xit) = 
(4/7r)/(f) sech^ ^(t). The dressing x which is not a CV 
internal mode can be observed as a modulation of the 
structure of the kink a. We proved that the existence 
of a directed energy current arises from the existence of 
the internal degree of freedom, T{t), combined with the 
dressing x{t). We fomid that the directed energy cur- 
rent vanished when F was set equal to Fq. When we set 
X = in Eq.(lO) for P the right hand side vanishes and 
we obtain 

P + /JP = 27r/i 

The infinite time average of this equation vanishes when 
we use the fact that the thermal average of the initial 
value of P vanishes. Thus there is no directed current 
in the SG unless the slope depends on t and the kink is 

dressed by x- 

We observe in the computation of {P{t)) that the 
kink sees the heat bath only through the damping term 
—f3{P{t)) i.e., {P{t)) does not see the fluctuations of the 
heat bath. The reason is that when we represent the bath 
as a generalized Fokker - Planck equation and calculate 
(P(i)) the damping term contributes —(3{P{t)) because 
the damping is represented by However the fluctua- 

tion term is proportional to a second derivative d^/dPdP 
and thus gives a vanishing contribution to {P{t)). Note 
a fluctuation such as (P^) or {P{t)P) would see both the 
damping term and the bath fluctuations. 

In conclusion, we have proven that the symmetry 
breaking that leads to a directed energy current in the 
ac driven SG is generated by the existence of the time 
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P{t) 





FIG. 3: The energy current P{t) and the slope of the kink r{t) for the parameters to = 0.25, /3 = 0.15, ei — 0.16, 62 — ei/\/2 
and 9 = 1.61 — n. Both curves show the effect of two driving frequencies. The time average of P{t) is nonzero and r{t) has 
a multiple frequency oscillation about {T{t)) greater than To with an amplitude change Ar/{r) ~ 20%. P{t) has the units of 
momentum and r(t) has the units of inverse length. 

P(t) 





FIG. 4: P{t) and r(t) for the parameters to = 0.3, (3 = 0.2, ei = 0.3, £2 = ei/V3 and 6 = 1.61 - vr. {P{t)} is nonzero and the 
slope r(t) has a multiple frequency oscillation about {T{t)) > Fq with an amplitude change AF/(F) ~ 50%. P{t) has the units 
of momentum and T{t) has the units of inverse length. 



dependence of the slope, T{t), and by the dressing, 
In Ref.^3' Salerno and Quintero showed that a double 
SG showed ratchet behaviour. In Ref . [l^ . Marchesoni 
obtained a directed kink transport by the sine/) poten- 
tial for £2 = 0. Costantini et al. observed ratchet 
behavior in ac driven assymetric kinks. 



where il{k) = (1 - 
Hnearized SG are 



A:^)^/^ and the eigenfunctions of the 



MO = (2^) 



ik — tanh^ 



(A.2) 



Acknowledgments 

We would like to thank Sergej Flach for useful discus- 
sions on symmetry breaking. 



APPENDIX 

We calculate the dressing x generated by the ac driven 
fi{t) in the linear approximation i.e., the lowest order 
in ei, €2 and by using the Green's function of the lin- 
earized SG equation in the presence of a soliton Eq.(C5) 
of Ref. [14]. The formal solution for x is: 

/oo / \ —1 poo 

dk(n{k)\ MO dOrkiO 
-OQ \ / <J —OO 

) 

dt' smn{k){t ~ t')e-^^*~*"> fi{t') (A.l) 



The integral is: 



) J — OO 

/OO 
dC'sinfc^'tanh^' (A.3) 
-OO 



The first integral is an irrelevant constant which we can 
neglect. Integrating the second integral by parts we ob- 
tain 



-{i/k) I cos kO sech^ £,'dC -{i/k)F{k) ( A.4) 
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FIG. 5: P(t) and r{t) for the parameters to = 0.1, (3 = 0.02, ei = e2 = 0.03, and 9 = tv. {P{t)) = as is required by 61 = tt. r{t) 
has a two frequency oscillation about (F) > Fo with only an amplitude change of AF/{r) ~ 4%. The relatively weak response 
of P{t) and r(t) is due to the smallness of the driver e. P{t) has the units of momentum and T{t) has the units of inverse 
length. 



where F{k) = 7rfc[2 sinh(7rfc/2)] ^ which decays rapidly 
with large k. The k integration in Eq.(A.l) is 



Rc(-i) 



dk e'^^[ik -tanh ^]k-^F{k) 
cos k^F{k) dk 
(7r/2)tanhW dfc sinfc^ [sinh(7rfc/2)] 



(A 



(2/7r)( sech^ ^ - tanh^ C 



— sech^ ^ 

TT 



2 



where we have dropped the irrelevant constant (2/7r) 
The time integral in Eq.(A.l) is f{t) where 



In this paper we only consider values of lo which are ap- 
preciabely less than one. While n{k) = (l + fc2)i/2 where 
in our units the lower band edge has the value one. Con- 
sequently with the w's in this paper as in Ref.[l] , there is 
essentially no radiationof SG phonons generated by the 
ac driver but only a dressing of the soliton that is lo- 
calized on the soliton. The presence of the SG phonons 
would not qualitatively alter the symmetry breaking but 
for the ac driver frequencies used in this paper the SG 
phonons would not be observable because they would oc- 
cur only in very high orders of perturbation theory. Since 
F{k) decreases rapidly with increasing k and ^ 1 we 
can treat f{t) as effectively independent of k and equal 



f{t) = (ei/2)cos^t^ ^2 ^/,..^2 + ^2 



f{t) 



dt'&mn{k){t~t')e-^^'~''^h{t') 



-f (£2/2) cos(2wi + 



1 - 2w 



(ei/2) cosut 



+ (£2/2) cos(2wi -I- 



/32 + {n{k) - Lof 

n{k) + UJ 
^ (3^ + {n{k) + uj)^ 
n{k) -2uj 



^2 + (1 - 2ujY 
l + 2uj 



13^ + in{k) - 2cj)2 
n{k) + 2lu 



^2 + {n{k) + 2cj)2 



(A.6) 



/32 + (1 + 2cj)2 

Finally we have 

X=(4A)/(t) sech^e 
with f{t) given by Eq.(A.7). 



(A.7) 
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